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ON A NEW METHOD OF COMPUTING THE ROOTS OF 
BESSEL'S FUNCTIONS 

By "William Marshall 

The object of this paper is to show the application of a certain method 
to the computation of the roots of BessePs functions. This method which 
was employed by me in my Doctor's Thesis * in the computation of the roots 
of the functions of the elliptic cylinder, and which is generally applicable in 
similar problems is essentially this : The coefficients of certain asymptotic ex- 
pansions which are necessary in the computation, are not found by arithmetic 
processes, but are determined from recurrence formulas which arise when 
certain differential equations are formally satisfied . The results are of course 
not new. The roots of the Bessel's functions of order zero as determined 
from the semiconvergent expansion were first given by Stokes, t The roots of 
the Bessel's functions of order n are given by McMahan t and are found in 
Gray and Mathews. § 

The ordinary Bessel's equation, namely 

^♦:^ + (> -a '•«-■ 

can by means of the substitution U= XKT n (x) be reduced to the form 

A solution of this in the neighborhood of x = oo is 

(3) U= C\P cos (a- x) + Q am (a- x)\, 

*Diss. ZUrich, 1909, American Journal of Mathematics, vol. 31, no. 4, Oct., 1909. 
t G. (i Stokes, Trans. Cambridge Phil. Soc, vol. 9, part I, or Math, and Phys. Papers, 
vol. 2, p. 329. 

J Annals of Mathematics vol. 9, Jan., 1895. 
§ Treatise on Bessel's Functions, p. 241. 
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where C and a are integration constants* and where P and Q represent the 
following asymptotic series : 



f» i (4n»-l)(4n*-9) 
1.2(8a:) 2 



(*) 



(4n 2 - 1) (4n 8 - 9) (4n» - 25) (4n 2 - 49) 
+ 1.2.3.4 (8x> 4 



4^-1 _ (4n 2 - 1) (4n 2 - 9) (4n 2 - 25) 
" ~ Sx 1.2.Zj8xj* + " ' 

If Q = -»/— and a = tt, and if Jf and i|r are determined so that 

V7r 4 

(4') M cos 1^ = P and ili sin i/r = Q, 

then from (3) 

(5) U= CM cos (a — x — -^r), 

and as 17 = xV„(as), 

and the roots of J n (x) = can be calculated from the relation 

(6) *+± w -z- + = -*Z± w( h=l t 9,...). 

M and yfr have been calculated approximately and with great labor from 

the relations M = \lP 2 + Q>, + = t*nr l (Q/P) respectively. They may be 

obtained much more simply by the aid of equation (2) as follows. 

U= CM cos (a - x — <r) wil1 satisfy (2) whatever values maybe given to C 

and a. Let C = 1, and a = 0. Then a Bolution of (2) is 

(7) U=M cos (yfr + x), 

and M and f are to be determined so that (2) is satisfied. From (7) we obtain 

(8) V = cos (^ + x)M' - MW + 1) sin (-f + *), 

* Gray and Mathews, 1. c, p. 87. 
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(9) U" = cos (>fr + x) M" - 2M' sin (f + x) (f + 1) 

- M cos (>fr + x) (f + l) 2 - M sin (-f + x)f". 
Substituting the values of U" and Cin (2) we obtain 

(10) M" cos O + x) - 2M' sin (f + x) (f +1)-M cos (-^ + x) (^' + l) 2 

1 5-=) = 0. 

Since the coefficients of sin (^ + x) and cos (^r + a;) must separately vanish, 
we obtain the two following differential equations for M and ^r, namely, 

(11) 2M'W + 1) + My = 0, 

(12) M" - mw + iy + m(i - ^^) = 0. 

Equation (11) integrates readily into 

(13) *' + 1 = i' 

where k is the constant of integration, which may be determined from the fact 

O PQ' — PO 
that ->fr = tan -1 -p» whence yfr' = — ^ -^- . Therefore h = 1 and we have 

(14) *' + 1 = ^ 

Substituting this value of i/r' in (12) we have a differential equation for M, 
namely 

(15) M»--± + M(l-^)=0. 

Now from (4') we see that the first term ofilf must be 1 ; we therefore assume 
as a solution of (15) 

(16) jf«i + ^ + £ + £ + ... f 

where the A'a are undetermined constants. From (16) we obtain 

( i„ «._!£4 + ±*i& + I3i + .... 



156 MARSHALL [July 

Also, by the Binomial Theorem, 

(™) W 1 + x* + tf + tf> + ' 

where 

a 2 = — 3.4,, 

( 1 9) «,= - 5^- 7^, _ 3A> 

_ — 6Ji 2 or 6 — 8A t a t — 10A^a 2 

a s — o^3. 8 , 



or in terms of the A's alone 

Q% = — OXl-^y 

a t = 6AI - SA iy 
(20) a 6 = - 10^41 + 1 2u4 4 J. 4 - 3 A 6 , 

a, e = 15 At- ZOAIA,, + 12 A t A 6 + %A\ - SA e , 



Substituting now in (15) the value of M from (16), M " from (17), and 
1/M Z from (18), we have 

(2n 3.2 A 2 5.4 A, 7.6 A 6 

^a 2 . -^4 . A e A s 
x 1 x* ar X s 

n* - i (« 2 - i) A )n» - j) A (n« -J) 

— -^i 6 — • • • 



8 8 
x 4 * 6 X s 



a2 a 4 a 6 cr 8 . . . _ o 
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Equating to zero the coefficients of like powers of x and making use of the 
relations (19) we have 

4A 2 = « 2 - i, 

<L4 4 =(n 2 -i-3.2)A--^p, 

(22) tA 6 = (»■_*_ 5.4)^ - 5M+ 3 Hfl3 , 

6u4 3 tf 6 + 8^4 4 a 4 + lOA 6 a 2 



4A = («» - J - 7.6)^ 6 - 



Formulas (22) then enable us to compute the coefficients in the expan- 
sion of M = y'P 2 + Q? without serious labor. These formulas are, in terms 
of the A' a, 

4A 2 = n* - 1, 
AA^ (nt-^Ai + GAl 

(23) 4A 6 = (n 2 -^A,- 10 Al + 12A 2 A it 

±A S = (n? - -^) A + 15 Al - ZOAlAt + 12 A^A* + 6A\, 

or in terms of n alone, _ 

A = ~ (4« 2 -l), 

A i= I 9 (4w*-l)(20n 2 -53), 

(24) A = -is ( 4m ' 2 - l)(240w 4 - 2488n 3 + 4477), 

A e = ~ (4n 2 -l)(12480)i 6 -322192n 4 +2023460w»- 3066403) 

it 
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or, if we set for shortness 4>< 2 = m, 

1 
^2 = ^r( m ~ I'). 

A= ^- (m -1) (5m -53), 

(25) J = -^ (m - l)(15m* - 622wi + 4447), 

A = -oiF (»» - 1)(195 m 3 - 20137?re 2 + 505865m - 3066403), 



From these coefficients in the expansion of M, we may readily compute 
those in the development of yfr = tan -1 $/P. We have from (14) 

(26) t'=j^-l- 
Then i/r' is of the form 

(27) +' = - 2 + fj [ + % + ■ . . . 
v r X 2 X* x e 

Integrating, and setting the constant of integration equal to zero, since 
tan ifr = QjP, we have 

(aa\ . _ ~ ^2 _ W_ _^_6 _ Cg _ 

K J T ~ x 3x 3 5x s 7a.' 7 



where the C"s are simply the coefficients in the expansion of 1/M 2 by the 
Binomial theorem, that is, 

20 i = -3A i C 2 -4A i , 
(29) 3C, = -4A a t\-5A l C,-6A 9t 

4C 8 = - 5A 2 C 6 - GAiCt - 7A 6 C 3 - 8A S , 
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METHOD OF COMPUTING ROOTS OF BESSEL'S FUNCTIONS 159 



or, in terms of the A's, 

a e = - 2A„ 

(30) C 6 = - ±A% + QA. 2 A t - 2A 6 , 

C a = 5 At - \2A\Ai + QA 2 A 6 + ?>A\ - 2A 6 , 

These coefficients in the expansion of yjr are used in the computation of 
the roots of J u (x) in the following manner : an asymptotic value of J n (x) is* 

(31) M x)=^[r C o S (^±±*-x) + G . in (*Lti «■-«)]. 

This may be written 

(32) J n (x) = ^A Jf cos (?^±1 „■-*_*), 

which vanishes when 

,,„, 2n+l 2ft - 1 

(33) — - — ir-x- + = ^ — ""» 

that is, when 

(34) x = 7 L(2n - \ + U-) - +, 

where k is any integer. 

7T 

If we put for shortness - (2» — 1 + 4i) = fi, 

then the values of x for which J„(x) vanishes are found by solving the 
equation 

(35) x = /8 - Vr, 
or 

(36) * = /3+^ + .^+ft+---- 

x ox ox 

Equation (36) is best solved by the method of successive approximations. 
The successive approximations are 

* Gray aud Mathews, Treatise on liesfel's Functions, p. -10. 
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(37) 



x 1 = fi lt 



C 2 



** = + j , 



*. = » + % + 



C A C 

— o 2 &Vj<i— - o 2 o 4 H — — 



P 



+ 



3 



£ 6 



* 5 = /3 + 



fv fy j 

W p 2 H« _ * /> /-y , arya 

C, 3 ~ ° 2 5 3 3 C * + 2C ' 



/8 + P 

7 



+ 



8 



^ 



Expressing these coefficients of /3 _1 , /8~ s • • • in terms of to by means 
of the relations (25) and (30) we have for the Ath root of J n (x), 



(38) 



x k = /3 



where 
and 



to - 1 (to - 1 ) (7to- 31) 
~Wp~ ~ 3-2 7 P 

(to - 1 ) (83m 2 - 98 2to + 3779) 
15. 2 10 ' /3« 

(to-1) (6949to 8 - 153855to 2 + 1585743to - 6277237) 
~ 105. 2 15 ' /3 7 - • • • 

/8«^(2n + 4*-l), 

to = 4n 2 . 



This result agrees with that given by McMahan ;* there is however a 
slight error in the formula as given by Gray and Mathews,t due to the drop- 
ping of a figure in the coefficient of yS — 7 . 



Pukdue University. 

OCTOBEK, 1009. 



* Annals of Mathematics, ser. 2, vol. 9, p. 25. 
fLoc. cit., p. 241. 



